Instructions on how to obtain Data_8C

1) Get el-ph coupling constants
 1-1) execute a scf calculation by pw.x in QE code with in.scf
      $MPIRUN -np $NUM_CORE $PATH/pw.x < in.scf
 1-2) execute a phonon calcualtion by ph.x in QE/PHONON code with in.ph
      $MPIRUN -np $NUM_CORE $PATH/ph.x < in.ph
 1-3) execute a el-ph calcualtion by ph.x in QE/PHONON code with in.elph
      $MPIRUN -np $NUM_CORE $PATH/ph.x < in.elph

2) Get Coulombic interaction Kernel (doesn't affect the results but is needed to run scdft.x)
 2-1) execute a scf calculation by pw.x in QE code with in.scf
      $MPIRUN -np $NUM_CORE $PATH/pw.x < in.scf
 2-2) execute a non-scf calculation by pw.x in QE code with in.twin
      $MPIRUN -np $NUM_CORE $PATH/pw.x < in.twin
 2-3) execute a Coulombic Kernel calculation by rpa.x in QE/sctk code with in.rpa
      $MPIRUN -np $NUM_CORE $PATH/rpa.x < in.rpa

3) Get eigenvalues used in scdft calculation
 3-1) execute a non-scf calculation by pw.x in QE code with in.nscf
      $MPIRUN -np $NUM_CORE $PATH/pw.x < in.nscf
 3-2) execute a pp calculation by get_eigval.x in QE/sctk code with in.eig
      $PATH/get_eigval.x < in.eig
 3-3) add a constant C to eigenvalues in pwscf.save/eigval.dat file (generated by get_eigval.x)
      to shift the chemical potential in scdft calculation
      
4) execute scdft calculation by scdft.x in QE/sctk code with in.scdft
   $MPIRUN -np $NUM_CORE $PATH/scdft.x < in.scdft

5) extract DOS and Z values from a standard output of scdft.x
6) repeat 3-3), 4), and 5) with different C until obtaining all the data
